
source("R/fct_Orthant.R")
source("R/Data_Orthant.R")
	
t<-proc.time()
reso<-order(A,a,m,M)
tOrd<-(proc.time()-t)


t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,2)
mi2<-proc.time()-t+tOrd
t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,3)
mi3<-proc.time()-t+tOrd
t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,4)
mi4<-proc.time()-t+tOrd

t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,5)
mi5<-proc.time()-t+tOrd
t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,6)
mi6<-proc.time()-t+tOrd





M3<-floor(mi2[[3]]/mi3[[3]]*M)
M4<-floor(mi2[[3]]/mi4[[3]]*M)
M5<-floor(mi2[[3]]/mi5[[3]]*M)
M6<-floor(mi2[[3]]/mi6[[3]]*M)

cat("\n")
cat("\n")
#source("R/test_Orthan_S.R")
#solst<-jst(A,a,m,M/10,1,1,10)
#
Mm<-50
#XX<-rtmvnorm(10000,mean=m,sigma=A,lower=a,upper=rep(Inf,p))
i=1
b<-matrix(0,nrow=Mm,ncol=5)
for(i in 1:Mm)
{
	b[i,1]<-mt(reso$S,reso$a,m,M,2)
	b[i,2]<-mt(reso$S,reso$a,m,M3,3)
	b[i,3]<-mt(reso$S,reso$a,m,M4,4)
	b[i,4]<-mt(reso$S,reso$a,m,M5,5)
	b[i,5]<-mt(reso$S,reso$a,m,M6,6)
	
	cat("loop ")
	cat(i)
	cat("\n")
}
z<-factor(c(rep("M2",Mm),rep("M3",Mm),rep("m4",Mm),rep("M5",Mm),rep("M6",Mm)))
boxplot(as.vector(exp(b[,1:5]))~z)
